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Abstract. In this manuscript we address the problem of deriving analytic expressions 
for calculating universal decoherence-induced errors in qubits undergoing arbitrary, 
unitary, time-dependent quantum-control protocols. For a qubit undergoing unitary 
decoherence the evolution of a qubit state in the presence of time- varying semiclassical 
fields may be treated geometrically. We show that the fidelity of an arbitrary control 
operation may then be expressed to arbitrary order in terms of experimentally relevant 
spectral characteristics of the noise and the control over all Cartesian directions and 
accounting for noise cross-correlations. We formulate control matrices in the time 
domain to capture the effects of piecewise-constant control, and convert them to 
generalized Fourier- domain filter functions. Such generalized filter functions may 
therefore be derived for complex temporally modulated control protocols, accounting 
for susceptibility to rotations of the qubit state vector in three dimensions. Taken 
together, this framework provides a computationally efficient means to calculate the 
effects of universal noise on arbitrary quantum control protocols without the need 
for time-consuming simulations of Bloch vector evolution. As a concrete example, 
we apply our method to treating the problem of dynamical decoupling incorporating 
realistic control pulses of arbitrary duration or form, including the replacement of 
simple TT-pulses with complex dynamically corrected gates. 



Arbitrary quantum control of qubits in the presence of universal noise 



2 



1. Introduction 

The importance of efficiently controUing the evolution of a quantum coherent system has 
given rise to the development of the field of Quantum Control Theory [U EJ [3l HJ [5l EJ [7j . 
Significant research attention has been devoted to this discipline due to the remarkable 
promise of quantum coherent technologies for future applications including quantum 
computation, quantum-enabled sensing, and quantum-enhanced metrology. 

A primary challenge being addressed by the research community has focused on 
understanding and mitigating the deleterious effects of decoherence on quantum systems. 
In these processes, uncontrolled interaction with the environment leads to randomization 
of a quantum system's state, effectively destroying its "quant umness." This phenomenon 
is especially important in quantum information settings where net error rates deep below 
fault-tolerance thresholds 10~^) are required in order to build scalable quantum 
computers [8]. 

A variety of methods have been developed to characterize the performance 
of quantum control protocols, with an eye towards estimating error rates due to 
decoherence in realistic, experimentally relevant noise environments which often have 
significant Non-Mar kovian characteristics. One of the most interesting from a practical 
perspective is the concept of spectral overlap formalized first by Kofman and Kurizki 
[9l [lOj. In this general approach, the net susceptibility of a given quantum control 
protocol to environmental noise is given by the overlap in frequency between the 
environmental noise power spectral density and the spectral characteristics of the 
modulation imparted by the control. 

These insights have been particularly important for the field of dynamical error 
suppression, which seeks to provide error robustness to quantum hardware at the 
physical level. These techniques address both implementation of memory or the identity 
(dynamical decoupling) [IT[ [121 1131 El 113 well as nontrivial quantum logic gates 
(dynamically corrected gates) [161 [IT] exploit interference trajectories to effectively 
time-reverse the accumulation of error. 

The concepts of the spectral overlap were expanded in the context of dynamical 
decoupling by Uhrig and Cywinski ^Hl 119! , in attempting to understand the efficacy of 
these protocols in suppressing decoherence due to Non-Mar kovian environments. It was 
shown that, in considering DES broadly, the relevant control protocol could be thought 
of as a noise filter^ and the efficiency of a particular control sequence deduced simply by 
an examination of the relevant filter transfer function [20j. This general approach has 
been repeatedly validated in experiments using a variety of technologies - from trapped 
atomic ions to semiconductor spin qubits [211 ESI ESI EH ES] . Remarkable agreement 
has been demonstrated between simple theoretical expressions for error derived from 
the so-called "filter-function" and experimentally measured noise power spectra - so 
much so that inverting this approach has allowed DES protocols to be used in noise 
spectroscopy [261 [271 [28] . 

Despite these capabilities, there remains a significant gap in our understanding of 
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how to efficiently calculate expected operational fidelities and error rates for complex 
quantum control protocols and for situations in which universal (i.e. multi-axis) time- 
dependent noise sources are present. For instance, relatively little is known about how to 
effectively account for the accumulation of error due to pulse non-idealities in dynamical 
decoupling sequences; the bang-bang limit is still widely assumed in analytic treatments. 

Generally, while it is known that spectral overlap techniques may be employed in 
order to evaluate error rates, it is not known how to efficiently produce analytic filter 
functions appropriate for such complex quantum control protocols. This is largely due 
to that fact that the analytic complexity of deriving such functions grows significantly as 
soon as non-commuting operators appear in the control Hamiltonian - as would be the 
case in dynamical decoupling with nonzero-duration control pulses or other nontrivial 
control operations. 

In this manuscript we address this challenge by developing a generalized method 
for evaluating the effects of universal, semiclassical decoherence on a quantum system 
undergoing an arbitrary unitary control protocol. By adapting the approach of Byrd 
and Lidar [29j, to the case of bounded controls in a time- varying noise environment, we 
demonstrate how the effect of universal decoherence on an effective spin- 1/2 qubit may 
be given to a geometric interpretation using a control matrix. In the presence of weak 
decohering noise, we show how all Cartesian contributions to the resulting entanglement 
fidelity may be calculated using noise power spectral densities and cross-power spectral 
densities along with the Fourier-space representations of the elements of the control 
matrix. Further, we show how it is possible to simply construct generalizations 
of the filter function for arbitrary control protocols by assuming piecewise-constant 
controls. As an example of the utility of this approach we study dynamical decoupling 
incorporating both realistic "primitive" ttx pulses and more complicated dynamically 
corrected gates. We are able to derive first order filter functions for the incorporation 
of arbitrary control pulses into a dynamical decoupling sequence and validate the 
performance of dynamically corrected gates in mitigating pulse errors in these sequences. 

The remainder of the manuscript is organized as follows. In Sec. [2] we provide 
the relevant background to the problem of decoherence in qubits. This is followed 
by Sec. |3] which introduces the theoretical framework for quantum control in noisy 
environments and a geometric interpretation of decoherence, leading to the derivation 
of analytic expressions for the fidelity in Sec. |4} Sec. [5] then provides a concrete method 
to evaluate the control matrix (and hence the generalized filter functions) for the case 
of piecewise-constant control in the presence of universal and pure-dephasing noise. As 
an example demonstrating the practical utility of our method, this approach is applied 
to the problem of finite-pulse effects in dynamical decoupling in Sec. [6} Our framework 
permits us to efficiently separate out contributions to the filter function from the pulse 
locations from effects arising from the pulse form in a dephasing environment. 
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2. Background 

A qubit is quantum system that exhibits two possible outcomes for the measurement of 
a particular physical observable: the z-component of the intrinsic angular momentum of 
a spin-1/2 particle being the archetypal example. Letting |0) and |1) denote the states 
in which one or other of these outcomes are returned with certainty, and assuming 
maximal knowledge of the system, an arbitrary qubit state may be written as a linear 
combination 

= cos(i9/2)|0) + e^^sin(i9/2)|l) (1) 

The angles (f) define a three-dimensional unit vector, the Bloch vector, that provides 
a convenient pictorial representation of the state [8j. Interpreting all possible states 
geometrically in this way, the two-dimensional qubit state space is mapped to a sphere 
of unit radius, called the Bloch sphere (Figure [T]). 

The temporal evolution of the qubit is described by the propagator U (t, 0), a unitary 
operator that transforms some initial state \ip{0)) to a final state \ip{t)) = U{t^ O)|'0(O)), 
for t > 0. The propagator satisfies the Schrodinger operator equation {h = 1) 

ijU{t,0) = Hit)Uit,0) (2) 

where H{t) is the Hamiltonian operator, which we assume to be traceless (i.e., 
Tr {H{t)) = 0), for all t. Any solution to § may be written in the form U (t, 0) = 
exp [—i9{t)n{t)a/2]^ where n{t) = {nx{t)^ny{t)^ Tiz{t)) is a real row vector of unit length, 
a = {crx^Gy^GzY ^ column vector comprised of the Pauli spin operators, and 6{t) is 
a real function of time. The action of U{t^O) on the state |'0(O)) can then be thought 
of as a rotation of the corresponding Bloch vector through an angle about an axis 
defined by n(t) . The set of all solutions to ([2]), in conjunction with the usual operator 
product, form a subgroup of the Lie group SU{2) [30j. 

A basic requirement for the realization of a variety of quantum-based technologies, 
in particular quantum information processing (QIP) devices, is the capacity to 
manipulate coherent quantum states of the form ([T]) with a high degree of precision 
in the presence of uncontrolled environmental interactions. For a quantum system 
even relatively weak interactions with the environment can have a dramatic effect. An 
initially coherent quantum system S and its environment E become entangled^ forming 
a single composite quantum system SE. From the perspective of an observer, with 
the capacity to measure only coherence is eflFectively lost [31j. If unaddressed, this 
process of decoherence can occur very rapidly, making the long term retention and 
accurate processing of information impossible. 

In this work, rather than employ the full quantum-mechanical approach described 
above, we model the effect of the environment through a fiuctuating classical noise field. 
Since the classical noise model and the purely quantum model of system-environment 
entanglement predict similarly negative impacts on the information processing capacity 
of a quantum system, we generically refer to both these processes as decoherence. The 
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semi-classical approach has the advantage of not requiring a detailed knowledge of the 
system-environment interaction Hamiltonian. Further, it has been shown repeatedly 
that such phenomenological models can accurately represent important aspects of 
experimental reality [22j. 

3. Quantum Control in Realistic Environments 

In this section we outline a mathematical treatment of a controlled quantum system 
affected by classical fields that serve to randomize the phase and probability amplitudes 
of a the quantum state, while preserving the vector norm (i.e. we omit leakage from the 
computational state space). 

We consider an ensemble (either spatial or temporal) of identically prepared qubits, 
with dynamics determined by the single-qubit Hamiltonian H{t) = Ho{t) + Hc{t). Here, 
the effect of the environment is modeled by the generalized noise Hamiltonian 



where /3{t) = {I3x{t)^ /3y{t)^ /3z{t)) is a three-element row vector, each component of 
which is a random process modeling classical noise in one of the three spatial dimensions. 
Control over the state of the qubit can be achieved through the application of an external 
field, represented by 



where the components of the vector h{t) = {hx{t)^hy{t)^hz{t)) describe the strength 
and direction of the control field as a function of time. We note that, for simplicity, 
we assume the all control is provided by Hc{t)^ with no contribution from the internal 
qubit Hamiltonian. 

3.1. Capturing the effects of time- dependent noise 

As a brief qualitative illustration of the effect of noise on the capacity to manipulate the 
state of a quantum system, consider the application of a simple control Hamiltonian 
Hc{t) = to an otherwise isolated qubit, where Q = tt/t for some positive 

real number r. Under these circumstances, the total Hamiltonian is independent of 
time and equation ^ is easily integrated to give [/(r, 0) = exp[—i7rax/2]. Thus, 
over a time interval [0, r] the control field generates a rotation of the qubit Bloch 
vector through an angle of tt radians (a 'vr-pulse'), about the x-axis (Figure [T^). If 
we allow for nonzero but constant 'noise' by letting /3{t) = (0,0,/?^), for example, then 
[/(r, 0) = exp[—i7r{ax + i^crz)/'2]^ where u = 2/3z/^. The propagator now describes a 
rotation through an angle 0' > tt, about an axis that is tilted away from the xy-plane 
(Figure [l|3). The noise therefore causes both a change in the angle of rotation and a 
shift of the rotational axis, so that the final qubit state will not be the one intended. 

In most cases, deriving such a simple closed- form expression for the propagator 
U{t) is not possible, as the Hamiltonian, in the presence of noise, will not generally 



Hoit) = f3{t)a 



(3) 



H,{t) = h{t)a- 



(4) 
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Figure 1. The Bloch sphere: showing a 7r-rotation about the x-axis, transforming 
an initial state |'0(O)) to a final state |'0(O)), for a) the ideal and b) the noise- affected 
control operations. The thick arrow indicates the final qubit state, the solid black 
arrow indicates the rotation axis, and the dotted path indicates the rotation angle. In 
b) the rotation axis and angle are both changed by the presence of the noise. Red 
extension of the path indicates overrotation in path 0'{t). See text. 



commute with itself at different times (from here on, we omit the first argument of 
a propagator when it is zero). Therefore, rather than attempt to solve ^ directly, 
we introduce the 'control propagator' t/c(t), defined as the solution to the noise-free 
Schrodinger equation idUc{t)/dt = Hc{t)Uc{t). We may then write U{t) = Uc{t)U{t)^ 
with the 'error propagator' U (t) capturing any deviation from the control evolution due 
to noise [32l |3^. Substituting U{t) Uc{t)U{t) into ([2]), and defining the 'toggling 
frame' Hamiltonian 

H^{t) = Ul{t)Ho{t)U,{t) (5) 

we find that U (t) satisfies the modified Schrodinger equation 

zfuit) = Ho{t)U{t) (6) 

In the absence of noise, executing a target operation Q over a time interval [0, r] 
simply requires one to engineer a control Hamiltonian such that Uc{t) = Q. Otherwise, 
the capacity to accurately 'steer' the qubit state via the control Hamiltonian Hc{t) is 
reduced by the presence of the noise Ho{t). In this case, we retain the requirement that 
Uc{t) = Q, so that U{t) = QU{t). 

3.2. Geometric treatment of decoherence: The error vector and control matrix 

The error propagator U (r), being an element of SU (2), may be interpreted as a rotation 
of the qubit Bloch vector through some angle 2a(r), about an axis a(r). Hence 

U{t) = exp [-^a(r)o-] (7) 
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where a(r) = a(r)a(r) is the real- valued 'error vector' with components a^(r), for 

1 /9 

i G {x^y^z}^ and magnitude |a(r)| = (a(r)a^(r)) 

That a geometrical interpretation of the action of a propagator is possible in terms 
of a rotation of the Bloch vector is a consequence of the homomorphism that exists 
between SU{2) and the group of three-dimensional rotation operators, S0{?>) [30]. In 
short, each propagator U and its physically equivalent negative — [/, both elements of 
SU (2), can be represented by a single rotation operator in S'0(3), in a way that preserves 
the group structure. 

We may use this geometric interpretation to understand how the effect of arbitrary 
control modifies the evolution of the Bloch vector in time, in the presence of 
environmental noise. From Eqns. we see that the error vector is determined 

by the toggling frame Hamiltonian HQ{t)^ which itself derives from the action of the 
control propagator Uc{t)^ an element of SU{2)^ on Ho{t). The homomorphism between 
SU{2) and S0{3) enables us to identify Uc{t) with a three-dimensional 'control matrix' 
R{t) G S'0(3), following [29j but now explicitly incorporating time-dependent processes, 

Hoit) = J2 muHt)a,U,{t) (8a) 

i=x,y,z 

/3i{t)Rij{t)aj (86) 



= f3{t)R{t)a (8c) 

Using the orthogonality of the Pauli operators with respect to the Hilbert-Schmidt inner 
product, the elements of the control matrix are 

R,,{t) ^ [R{t)],, = Tr {uKt)aMt)a,)/2 (9) 

for ij G {x,y,z}. 



4. Calculating the Operational Fidelity 

Our goal is thus eflFectively to characterize the argument of Eq. [7| using the formalism 
introduced above. In the specific case of dynamical error suppression, we go further 
and attempt to minimize this term. Accomplishing these goals will involve a series of 
analytic advances in how we calculate operational fidelity and capture the action of the 
control Hamiltonian. In the following sections we will describe our method in detail. 



4>1> Approximating the error vector using the Magnus expansion 

We seek an expression for the error vector a(r) by putting the solution to ^ in the 
simple exponential form ([7| . Due to the non-commuting nature of the operators aflFecting 
the quantum state during its evolution, this is not a straightforward task. However, at 
the cost of expressing the exponent as an infinite series, we can achieve the desired 
outcome using the Magnus expansion ^3ll [35l |36] . Specifically, given equation ^ , we 
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may write U{t) = exp [— i$(r)], where 

oo 

$(r) = ^$^(r) 



and 



(r)= r dtHoit) 
Jo 



$1 

$2(r) = - 
$3(r) = - 



dti 
dti 



dto 



dto 



Ho{ti)^Ho{t2) 

t2 



dt:^ 



^o(ti), [Ho{t,),Ho{ts) 



(10a) 

(106) 
(10c) 

(lOd) 



with higher order terms involving increasingly complicated multiple integrals of nested 
commutators. Here the index fi indicates the order of the Magnus expansion. 

We now exploit this general approach and the geometric interpretation of control to 
find expressions for the error vector. Returning to the control matrix introduced above, 
we compactly express Eq. ( [8c| ) as 

Hoit) = J2 A(^)n(t)<T (11) 

i=x,y^z 

where we have written R{t) = {rx{t)^ry{t)^rz{t)Y^ with r^(t) = {Rix{t) ^ Riy{t) ^ Riz{t)) 
a row vector comprised of the i-th row of the control matrix. Substituting ( [TT| ) 
into the above expressions for the Magnus expansion, and using the vector identity 
[ucr, V(t] = 2i(u x v)cr, for u, v G M^, we are able to write U{t) = exp [— ^ ci^(^)cr], 
with the error vector in the form of an infinite series. The first three terms of this 
expansion are 



ai(r)= [ dt(3,{t)T,{t) 



a2(r) 
a3(r) 



E 

2 

3 



dti 



tl 



dti 



(12a) 
(126) 

dt3/3i{ti)/3j{t2)/3k{t3)rijk{ti,t2,t3) (12c) 



dt2/3i{ti)/3j{t2)rijiti,t2) 
dto 



i,j,k=x^y^z ' 

where we've introduced the vectors r^j(ti,t2)) = [^i{ti) x ^jih)] and fijkih^h^ts) = 
[r^(ti) X [rj{t2) X Ykits)]] + [vkits) x [r^-(t2) x r^(ti)]]. Generally, the n-th order 
term is an n-fold integral over the sum of all possible products of the form 
Ai(ti)A2(*2)...An(^n)nii2...in(^i^^2, - - - , tn) , where the vector ri^i^...in{ti,t2, ...,tn) is a sum 
of multiple vector cross products of rows of the control matrix R{t) , evaluated at times 

tl, •••5 tn- 
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4^2. Expressing the fidelity 



One way of measuring how well a given propagator Vi approximates a target operation 
V2 is to take the Hilbert-Schmidt inner product ^1) = | Tr (V2^1^i) which, by analogy 
with the usual inner product of state vectors, effectively measures the 'overlap' between 
the two operators. Taking the square modulus of the inner product of U{t) = QU{t) 
and (5, gives J-'(r) = || Tr(C/(r))p, an expression equivalent to the so-called single-qubit 
entanglement fidelity [371 EH] • This then serves to quantify the accuracy of a control 
operation, subject to a particular realization of the noise vector /3{t). Experimentally, 
however, it is usually only the ensemble average of the entanglement fidelity 

J^Ut)^\{\Tv{U{t))\') (13) 

taken over all realizations of the noise Hamiltonian ([3]), that is measured, and we 
therefore refer to (13) as determining the 'fidelity' of a particular control operation. 

With the error propagator written in the form ([7|, and suppressing the explicit 
r-dependence, one finds 



^a. = ^[(cos(2a)) + l] 



(14) 



so that, for each qubit in the ensemble, the fidelity is determined by the magnitude of 
the angle through which the noise causes the control operation to be rotated away from 
the target operation Q. Expanding the cosine term in a Taylor series we obtain 

x-ir 



nil 



1 



1 + 



CO 

m=0 



,2m\ 



(2m)! 



(15) 



We proceed by evaluating the first nontrivial term (m = 1) in the expansion (the terms 
(a^^), for m > 1 are obtainable similarly). Recalling that = aa^, and using the 
Magnus expansion of the error vector ( |12a| ), we may now write 

{a') = 5](a,an = [(a?) + (al) + ... + 2 ((a^a^ + (aiaj) + (a^a^) + ...)] (16) 

where = aiaf , and /x, u are indices over the order of the Magnus expansion. Finally, 
we arrive at the following series expansion for the fidelity 



T =1 

nil J- 



(al) + 2(aiaJ) 



with al = 3-2^2 



al) - 2(aia^) ■ 

(aiaf)^. The second term in (17) 

/ dt2 / rfti(A(ti)/3,(t2))n(ti)rJ(t2) 



(17) 
(18 a) 

hj,k=x,y,z ' 

contains all two-point noise cross-correlation functions (A(^i)/3j(^2)), for hj ^ V-, ^}-> 
while the third contains all those evaluated at three time points, and the terms in square 
brackets all involve four-point correlation functions (this is determined by the sum of 



V dt2 dti{/3iiti)/3jit2))Rikih)Rjkit2) 

u-^..Jo Jo 
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subscript indices in a term as they indicate the expansion-order of the error vector in 



Eq. 12a 12c). For noise that is weak, the general trend wiU be for terms involving 
higher-order correlation functions to have less of an effect on the fidelity. Explicit forms 
for these terms, in the case that noise is purely dephasing (i.e. Ho{t) = /3z{t)az)^ may 
be found in [39j . 

4.3. Moving to the Fourier domain 

The form of the expressions contributing to the fidelity above may be calculated 
explicitly, but rely on time-domain correlation and cross-correlation functions and 
convolution with the multidimensional control matrix. Relying on the properties of the 
Fourier transform it is experimentally convenient to write these terms in the frequency 
domain. Generally, we can define the Fourier transform Si^,,,i^{uJi^ of an n-point 

cross-correlation function via 



(2vr) 

The fidelity ([l7|) can then be rewritten as 



n=2 \ i\...irL ) 



T =1 

where 7^ii..i^(cL;i, ...,cl;^) is determined solely by the control matrix. 



For noise that is sufficiently weak, the n = 2 term in (20) will be the dominant 
contributor to fidelity loss in the short term, i.e. J^av ^ 1 — {af). If the noise is also 
wide sense stationary then the two-point cross-correlation functions depend only on the 
time difference ^2 — ^i, so that 

(A(ti)/3,(t2)) = — / dujS,,{uj)e^-^'^-'^^ (21) 

where Sij{uj) is the cross-power spectral density between the random variables /5^(t) and 
(3j (t) . We then have 



2 ^ 



[cv) (22) 



ij,k=x,y,z 



where 



Rij{uj) = -iuj dtRij{t)e''^^ (23) 
Jo 

are the elements of the control matrix in the frequency domain (again we suppress the 



explicit dependence on the total time r). The form of (22) is reminiscent of the filter 
function introduced in previous literature on dynamical error suppression, to which 
we will return later. Thus, in the weak noise regime, it is possible to understand the 
fidelity of a qubit's unitary evolution in terms of spectral properties of the noise and 
the control [9l [I8] . In the next section we will show how the control matrix may be 
evaluated for the paradigmatic case of a piecewise constant control sequence. In doing 
so, we explicitly demonstrate the connection between our generalized approach and the 
well known filter-function formalism. 
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Qn 



Qn-1 



II 



Qn-: 



II 



II 



n-2 



// 



Q, 



t2 



Pi 



to 



Time 



Figure 2. Explanatory diagram for a piecewise-defined control sequence. During each 
interval a control operation Pi is executed. The operators Qi are the control 

propagator evaluated at the start of the l-th pulse, for / = 1, 2, .., n. 



5. Evaluating the control matrix 

In general, for arbitrary continuous-time modulation of Hc{t)^ the control matrix can be 
evaluated only approximately. However, using the simplifying assumption of piece wise- 
constant control it becomes possible to find an exact analytic form for the control matrix. 
This model can be applied to dynamical decoupling sequences (vr-pulses separated by 
free-evolution periods, ultimately implementing the identity operation /), as well as a 
variety of nontrivial control operations such as piecewise-defined dynamically corrected 
gates and composite pulses . 



5.1. Universal noise 

For a sequence of n consecutive unitary control operations Pi, P^, we divide the total 
sequence time r into n corresponding time bins such that Pi = Uc{ti^ti-i)^ for 

/ G {1,2, ...,n}, with tn = r and to = (Figure [2]). Letting Pq = /, and defining the 
cumulative operators Qi = P/P^_i...Po and Q = Q^^ equation ^ becomes 

= -5]<^^'H^)Tr (gj_if/l(t,t;-i)aif/e(t,t;_i)g;_ia,-) (24) 

1 = 1 

where the function G^^\t) = B[t — — Q\t — ti] has unit value within the /-th time bin 
and is zero elsewhere. Equation ([24]) can be rewritten in terms of the control matrices 
R^^{t — for each of the n operations Pi in the sequence, defined relative to the 
pulses start times by letting 

Ul{t,ti_r)cj,U,{t,ti_r) = (25) 

k=x,y,z 

for i G {x^y^z}. Again, the components of R^^{t — are found using the Hilbert- 
Schmidt inner product 

R^{t-ti_,) = = \Tv{U^^{t,ti_,)aiUe{t,ti_i)a,) (26) 



for i^j G {x^y^z}. Substituting ([25]) into (24), and using the linearity of the trace 
operation, we find that we can write the control matrix for the pulse sequence in the 
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compact form 

n 

R{t) = ^G^^\t)R''^{t - ti_i)A^^^ (27) 

where the matrix A^^^ has components 

A« ^ ^ Tr (QUa^Qi-iaj) (28) 
for i^j G {x, z}. Using ([23]) to switch to the frequency domain, we have 



R{uj) = ^e^"^^-^12^^(cj)A(^) (29a) 

1=1 

with 

R^^{uj) = -iuj ' ' dte'^^'R^'it) (296) 



JO 

With the above expressions and a knowledge of the statistical properties of the noise, 
using the results of Sec.|4]we can find an approximate value for the fidelity of an arbitrary 
piecewise-defined control sequence in a universal weak noise environment. Despite the 
proliferation of terms in these analytic expressions, this approach is straightforward to 
implement by hand or through numerics. Ultimately the computational efficiency of 
calculating the control matrix and spectral overlap integrals far exceeds that of brute- 
force numerical simulations of the evolution of the Bloch vector [39j. 

5.2. Pure Dephasing 

In many of the physical systems with the potential for use in nascent quantum 
information technologies, the characteristic single-qubit dephasing time T2 is 
considerably smaller than the relaxation time Ti, and governed by different physical 
processes. It is important, therefore, that we understand the effects of purely dephasing 
noise on single-qubit control. To this end, we consider a model system in which the 
noise vector /3{t) has only a z-component and the noise Hamiltonian ^ reduces to 

Ho{t) = (3,{t)a, (30) 

Only the third row of the control matrix, r^(t) = {Rzx{t) ^ Rzy{t) ^ Rzz{t)) ) is required 
to calculate the fidelity. So, for the sake of notational simplicity, we drop the z 
subscript from this point on and refer to r{t) = r^(t) simply as the 'control vector', 
with components r^(t) = Rzi{t)^ for i G {x, ^, z}. 



With the dephasing Hamiltonian (30), the lowest-order contribution (22) to the 



fidelity (17) simplifies to 



2^M^iH (31) 



where S{uj) = Szz{^)-> and the first order 'filter function' Fi{uj) is simply the square 
magnitude of the frequency domain control vector 

Fi(a;) = |r(a;)r= J] \rM\' (32) 



i=x,y,z 
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If we further assume that the noise is Gaussian, with a mean value of zero, then 
aU correlation functions for which n is odd are equal to zero. The 

Gaussian moment theorem also allows us to write all remaining higher-order correlations 
in terms of only the simplest two-point correlation function. In the frequency domain, 
the statistical properties of the noise are then captured entirely by the power spectral 
density S{(ju)^ and the fidelity (17) can be written as the infinite sum [39] 



^ n: 



1 - 



1 r doo 

^ / —S{u})Fi{uj) 



+ 



Here, the effect of the control on the qubit fidelity is described by the functions Fp, for 
p = 1,2, which depend only on the control vector r. The complexity of the filter 
functions increases rapidly with increasing order. However, Fi{(ju) can be quite simple 
in form and reasonably straightforward to calculate. In the weak noise regime we ignore 
the higher order terms in (33) and write 

^a.(r)~^{l + exp[-x(r)]} (34a) 

where the rate of fidelity decay is given by 

1 Hi J 

X{r)^{al) = - -S{uj)F,{u) (346) 

Thus, the decay rate is determined by a simple overlap integral between the power 
spectrum of the noise and a "filter function" representing the control in the Fourier 
domain. This result is similar to the familiar expression for the coherence of a dynamical 
decoupling sequence in the case of unbounded controls [18j (see also Sec. |6]). The 
difference here being that the filter function Fi captures only the lowest order nontrivial 
effect of the control on the fidelity in a dephasing environment, rather than providing 
an exact solution in the bang-bang limit. 

For a piecewise-defined sequence, Eq. (24) reduces to the following expression for 
the components of the control vector 



ri{t) 



1 ^ 

- 5] G^'\t) Tr (g;_it/](t, ti_,)aMt. 
1=1 



Equations (29a) and i\29bh become 



e^'^*'-ir^'(a;)A(') 



Pi I. .\ — 



-lUJ 



tl-tl- 



dte 



it) 



respectively, where the components of the l-th pulse control vector are 

1 



rf (t 



ti-i) = - Tr {Ul{t,ti.i)aMt,ti.i)aj) 



(35) 



(36 a) 



(366) 



(37) 
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for j G {x,y,z}. 

If we assume that during each of the n time bins comprising the sequence the 
apphed control effects a steady rotation of the qubit Bloch vector through an angle 6^/, 
about an axis fii = cos {(/)i)x + sin (0^)^ in the x^-plane, then for t G t/] 



Uc{t,ti_ 



exp 



-Qi {t - ti_i) a^f^i 



(38) 



where = 9i/{ti — t^-i) is the rotation rate and a^^^ = cos {(j)i)ax + sin {(l>i)cry. From (37) 



and (366) we can show that the frequency domain control vector for the l-th operation 



has components 



.Pi 



CaJ 



{S^j [i^igi{u)) - u)fi{u))] + [^ifi{u)) - ic^gi{c^)] Tr {a^^a^aj)} 



(39) 



where //(cj) = e^^(^^-^^ 



i)cos(^^) - 1 and gi{uj) = e^^(^^-^^- 



)sin((9/). Although ^ has 
been derived assuming rotations take place about axes in the x^-plane only, periods 
of free evolution (no rotation) and z-axis rotations are easily accounted for by setting 
9i = and f]/ = in (39), whenever one of these operations occurs in the sequence. 



This is possible simply because, for dephasing noise, both of these operations commute 
with the noise operator. 

With these tools, the process for finding the operational fidelity of an arbitrary 
piecewise-constant control sequence in the presence of dephasing noise may be 
summarized as follows: 



(i) Calculate the matrix given by Eq. (28) for each step in the chain of control 
operations 



(ii) Substitute the above result as well as Eq. (39) into Eq. (36a) 



(iii) Repeat for all time steps / G [1, n] to find the control vector in the Fourier domain 



(iv) Find the first-order filter function from Eq. (32) 



(v) Use Eqs. (34a) and (346) to calculate approximate expression for the fidelity valid 



for weak dephasing noise. 

In the next section we employ this approach to treat the important problem of 
studying pulse imperfections on the performance of dynamical decoupling sequences 
used in quantum memory. 



6. Example: Pulse effects in dynamical decoupling 

The aim of dynamical error suppression protocols, in the language described above, 
is to perform a desired control operation Q without significant error by constructing 
Hc{t) such that it modifies the action of the noise Hamiltonian to affect U{t) /, the 
identity operation. In particular, a large body of literature has been devoted to the 
study of dynamical decoupling (DD) sequences. These are control protocols aimed at 
implementing the identity operator {Q = I) and thus preserving the initial state of the 
qubit; decoherence is effectively time-reversed by the action of the control sequence. 
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Figure 3. Dynamical decoupling sequence for n identical pulses, each of width rp, 
centered at times ^^r, for / = 1, 2, .., n 



In the limit of instantaneous control pulses and in the particular case of a semi- 
classical dephasing environment (or equivalently a spin-boson model), it has been shown 
that it is possible to write the qubit coherence at the end of an arbitrary DD sequence, 
given preparation of the system in a superposition of eigenstates of a^, as Vl^ = e"^^'^^ 
[42[ E^- this expression, the effect of both the noisy dephasing environment and the 
control modulation are captured in the overlap integral 

Xir) = - [ ^S{u)F{u)du. (40) 



The terms in this integral include the statistical properties of the noise in the Fourier 
domain through the noise power spectrum S{(jj)^ and the action of the control sequence 



through the filter function, F{uj). This expression is formally identical to equation (346) 
above, with the exception that, because the control pulses are assumed to be unbounded 
in strength with vanishing duration (the so-called "bang-bang" approximation), the filter 



function in (40) is exact. The bang-bang approximation is, however, unphysical and has 
been shown in experiment to neglect important contributions to the net error arising 
from the nx control operations with nonzero width r^, that constitute the sequence (see 
Fig.g. 

In this section we will apply the techniques we have developed to evaluate the 
fidelity of dynamical decoupling sequences when considering physical constraints - such 
as the incorporation of realistic control pulses. In order to understand how well a 
given DD sequence suppresses error we must efficiently evaluate the filter functions, 
accounting for the various errors that arise when performing realistic control pulses in 
a dephasing environment. We adapt the results presented in section ^ for an arbitrary 
control sequence to the specific case of dynamical decoupling with bounded pulses of 
finite width, and show how the form of the pulses may be analytically separated from 
the applied pulse sequence. 

6.1. Separating pulses from the pulse sequence 

DD sequences can differ in both the spacing and form of the pulses in time. To find 
optimal error-suppressing sequences for a given environment or physical system, either 
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of these variables may be changed independently of the other, and these differences must 
be easily separated in an analytical framework. We are therefore motivated to express 
the sequence control vector r(cL;) in terms that can be separated into those dependent 
only on pulse location {5i) and those dependent on pulse type (expressed in terms of the 
pulse control vector r^(cL;)). 

As in section ([5]), we start with a piecewise constant sequence of n pulses, however, 
here we exploit the simplicity of the actual control Hamiltonian to streamline the 
process. Explicitly, we assume that the pulses are identical, each of width Tp and 
executing a net vr-rotation of the qubit Bloch vector about the x axis, i.e.. Pi = a^^ 
for / G {1,2, In addition, both before and after each pulse we include a period 

of free evolution (i.e., no control operation), so that immediately following the /-th 
pulse, we have Qi = PiIPi_iI ...IPq. In effect, we have two interleaved subsequences: 
a subsequence of n disjoint control pulses and a subsequence of n + 1 disjoint free 
evolution periods. Within the total sequence time r (including both control pulses and 
free evolution periods), the center of the /-th pulse occurs at a time 5it^ where < Si < 1 
and 5n^i = 1 (see Figure |3]). 

Due to the cumulative nature of the control propagator, the pulse and free evolution 
subsequences are not independent. In particular, the operator Qi results in change in 
sign of both the pulse and free evolution control vectors at the conclusion of each pulse. 



From (36a) and (39) (with Qi = 0)^ the control vector for a sequence of free evolution 
periods extending over the disjoint time intervals [0, 5ir — rp/2],[5ir + Tp/2^S2r — 
rp/2], [5nT + rp/2,r], with a change of sign during each successive interval, reduces 
to a single component 

n 

r,{uj) = 1 - e^^^ + 2cos {ujTp/2) J^(-l)^e^^^^^ (41) 

1=1 

The same result was derived previously by assuming finite-width pulses during which 
noise was ignored f22l [2T]. In the bang-bang limit of infinitely narrow pulses (r^ 0) 
we have the standard result (for even n) 

n 

r{uj) = 1 - e^^^ + 2 ^(-l)^e^^^^^ (42) 
1=1 

with F{uj) = \r,{uj)\^. 

Incorporating noise during pulses, arbitrary pulse forms, and arbitrary DD 
sequences has heretofore been impossible by simple analytic methods. By taking account 
of the terms dependent on r^, the accumulation of error during the control pulses 
and the effects of pulse shape on this error can be calculated for arbitrary decoupling 
sequences. Using a method similar to that above, we can derive the control vector for the 



subsequence of disjoint pulses - with arbitrary form - alone. Adding this to (41) gives 
us the control vector for a total decoupling sequence, which now has the two nonzero 
components 

n 

r,{uj) = 1 - e''^^ + [2 cos {u}Tp/2) - e'^'^^p/Vf (w)] ^(-l)'e*'^'^'^ (43a) 



Arbitrary quantum control of qubits in the presence of universal noise 



17 



ry{uj) = -e-^""-/Vf (cj) 5^(-l)^e^"^^" (436) 

These expressions capture modification of dephasing dynamics due to the duration and 
form of the pulses as weU as depolarization effects occurring due to pulse-errors derived 
from a fiuctuating detuning from resonance (i.e. the impact of dephasing noise). 

We will move forward using these expressions in order to treat experimentally 
relevant cases of DD incorporating various pulse forms including error-suppressing 
dynamically corrected gates. 

6.2. Pulse Forms: Standard ttx- Pulses and Dynamically Corrected Gates 

For a DD sequence we are interested in physical control operations that implement 
a logical NOT gate. The simplest approach is to employ "primitive" vrx-pulses, 
implemented by applying the ideal control Hamiltonian He = over a time interval 

[0,r7T-], where Q = n/r^. This induces a simple rotation of the Bloch vector through 
an angle of tt radians about the x-axis. With this geometrical picture it is possible to 
understand the accumulation of error during a pulse using the method we have described 
in detail above. In the presence of time- varying dephasing noise, the bare gate is modified 
by an error term that, to first order, results in a small additional rotation about an axis 



that has both y and z-components. From equations (36a) and (39) we find that the two 
corresponding components of control vector are 

rf = rf (uj) = (e-- + l) (44a) 

< = rf-)H = -^(e-- + l) (446) 

It is interesting to note the correspondence between the prefactors here and those arising 
in a master-equation treatment of a driven quantum system in the presence of dissipation 



Moving beyond the standard yr^-pulse of finite duration, we consider a gate that 
has been designed to provide robustness in the presence of random gate errors - the 
dynamically corrected gate. In previous work we demonstrated how such gates provide 
enhanced resistance to error in the presence of time- varying dephasing noise [39]. Here 
we include expressions for a dynamically corrected logical NOT gate consisting of three 
successive ttx pulses, the second of which is executed at a rotation rate half that 
of the other two ^44j- The control Hamiltonian for this sequence is 

0<t <T^ 

^(Jj^ r^<t< 3r, (45) 
naj2 3T^<t< 4r, 
In the case of dephasing noise, the control vector again has two components: 



^2 _Q2 ^2 _ (^/2)2 



(46 a) 
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Pi{uj) P2(^)/2 



(466) 



where pi{uj) = e^^^^^ + e^^^^^ + e^^^^ + 1 and p2{uj) = e^'^^" + e'^^". 
^.5. Interpreting the filter function 

The inset to Figure shows the filter functions, Fi{uj) = \ry{uj)\'^ + |rf (cc;)p, for 



the primitive and dynamicaUy corrected logical NOT operations described above, as a 
function of the angular frequency in units of r~^. The filter function grows with cj, 
capturing the low-frequency "filtering" performance of DD sequences. Beyond a certain 
frequency Fi{(jj) 1, above which noise is passed unimpeded. This corresponds to 
the physical observation that fiuctuations fast relative to the shortest interpulse period 
cannot be effectively suppressed by dynamical decoupling. 

The filter function's behavior near zero frequency captures the leading-order error 
susceptibility. We may express Fi{uj) oc (cjrT^)^^^^-^^ as (jUt^j; 0, then a is the order of 
error suppression for the sequence [45]. Using these expressions, an examination of Fig.ji] 
reveals that the relative performance of the primitive and dynamically corrected gates is 
captured in the slope of the filter function as plotted on a log-log scale ^20^ ; higher slope 
corresponds to a higher order of error suppression (for noise near zero frequency) due 
to the multiplicative nature of the filter function in ( |346 ). In a filter-design framework 



we may refer to this slope as a filter roUoff in the stopband, using the relation that the 
roUoff is 6{a + 1) dB/octave. At low frequencies the primitive vrx-pulse filter function 
is approximately quadratic, i.e. F{(jj) oc {ujTt^Y^ resulting in a = - the same as for free 
evolution. This is not a surprising result, since there is no error suppression built into 
the primitive gate. By contrast the DCG filter function scales as F{(jj) oc (c^Tt,-)^, giving 
a = 1, as expected for a construction designed to suppress error to first order. 

6.4' The full filter functions: DD -h DCG 

We now combine the control vectors derived for the pulses themselves with those 
obtained for the generic DD sequences, considering three paradigmatic cases: (1) bang- 
bang decoupling, (2) standard pulses with nonzero r^, (3) DCGs. This is accomplished 



by calculating and inserting the relevant Cartesian components into Eqs. 43a and 43 & 
for the full DD sequence. This study reveals how it is possible to evaluate complex 
perturbations to quantum control protocols in a straightforward manner using our 
method. 

We compare the error-suppressing performance of CP and UDD sequences as 
illustrative examples, accounting for each of the cases outlined above. The CP sequence, 
first proposed by Carr and Purcell in the context of NMR [46j, is a straightforward 
extension of the original Hahn spin echo [47j. In an n-pulse CP sequence, the /-th pulse 
has the fractional location Si = {I — 1/2) /n. The bang-bang limit CP filter function 
for n = 6 is shown as a black line in Figure lik) for n = 6. It exhibits an order of 
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error suppression of a = 2, as do all CP sequences with n > 2, giving a roll-off of 18 
dB/octave. This is a limiting case in which pulse effects are entirely ignored. 

If we now take into account the width of the pulses by assuming that they take the 
form of primitive NOT gates of width r^, as described above, then the components of 
the pulse control vector r^{(jj) are given by equations (44a) and (446). The resulting 
filter functions, indicated by the blue lines in Figure |4^) show a reduced order of 
error suppression a = 2 ^ 1, with a corresponding reduction of the filter roUoff to 
12dB/octave. As the pulse width grows as a fraction of the entire sequence, the relative 
importance of the pulses grows, as expected. This is manifested as an expansion of the 
range over which the modification to the filter-function due to pulse effects dominates 
the bang-bang filter function. 

Despite the reduction in the order of error suppression, there is still noise 
cancelation, regardless of the pulse width, since the sequence is in essence an Eulerian 
dynamical decoupling (EDD) sequence [48j. Replacing the primitive pulses with 
dynamically corrected composite gates and the components of the pulse control vector 
with (46a) and ( |466| ), we almost completely restore the original ISdB/octave roll-off for 
the CP sequences (Figure [3^). These filter functions lie approximately beneath the bang- 
bang filter function in this figure. This observation confirms that the DCG provides 
one order of error cancellation during the applied pulses. In essence, the deleterious 
effect of error accumulation during nonzero-duration control pulses can be mitigated by 
appropriate choice of a DCG. 

Studying UDD is more instructive because the n-pulse UDD sequence, with relative 




Figure 4. First order filter functions as a function of dimensionless frequency. The 
inset of a) compares filter functions for primitive (blue) and dynamically corrected 
(red) NOT gates while the main plot shows filter functions for a 6 pulse CP sequence. 
The results for primitive NOT pulses are shown in blue, while those for dynamically 
corrected pulses are shown in red. Plot b) shows the same set of results for the 6 pulse 
UDD sequence. Black lines represent filter functions for sequences in the bang-bang 
limit 
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pulse locations given by 5i = sin^ [7r//(2n + 2)], ensures that (in the bang bang limit) the 
first n-derivatives of rz{uj) vanish at = [I8l|42], giving an order of error suppression 
that increases linearly with n. The effects of pulse width and shape are therefore 
relatively more dramatic due to the high order of error suppression for, e.g. n = 6, as 
shown in Fig. |4]3). For sufficiently long pulses we see that the general benefits derived 
from the use of a UDD sequence can largely be destroyed as the error-susceptibility of 
the pulses dominates the suppression provided by the pulse timing. Again, we observe 
that the reduction in the order of error suppression arising from addition of nonzero- 
duration TTx-pulses can be partially offset by use of a DCG. However, the fact that 
by design, the DCG used here only provides first-order error suppression, reduces our 
ability to recover the original order of error suppression provided by a UDD sequence 
with large n. 

Nonetheless, these results indicate that judicious choice of DCG or other 
compensating pulse protocols within DD sequences provides a path to mitigate the 
effects of pulse errors in DD sequences. This is especially important in long-storage 
settings where large pulse numbers may be employed in order to effectively suppress 
dephasing errors [49j. These observations represent an important validation of our 
method, and show that it provides a straightforward approach to account for a variety 
of pulse-duration, modulation, and shape effects. 

7. Conclusion 

Understanding, predicting, and mitigating decoherence remain significant challenges 
for the development of quantum technologies. In particular, in order to effectively 
develop novel quantum control protocols that suppress error or evaluate the performance 
of a complex quantum error correction procedure, we require methods permitting a 
system designer to estimate error probabilities using real inputs about the environment. 
This capability is vital as the quantum science community moves from proof-of- 
principle concepts towards engineering of real quantum coherent technologies where 
rigid performance bounds are required. 

In this paper, we have presented a novel method for determining the fidelity of 
arbitrary quantum control operations in a universal time-dependent noise environment. 
Building on past work employing spectral overlap functions, our method provides a 
straightforward, easily automated approach to producing analytic functions capturing 
the effect of the control in the Fourier domain - generalized filter functions. Using these 
we are able to evaluate control fidelities to arbitrary order, for arbitrary time-dependent 
control operations in terms of experimentally relevant characteristics of the noise. While 
the error model we study is constrained to ignore, e.g., leakage errors from the qubit 
subspace or system-bath entanglement, the selected semiclassical model of universal, 
Non-Mar kovi an, time- varying noise is far more realistic than earlier approximations 
based on uncorrelated errors. 

We have exploited the capabilities provided by this technique to address the 
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key challenge of accounting for nonidealities in control operations used to implement 
error-suppressing dynamical decoupling sequences. By generating filter functions 
incorporating pulse errors and even the effect of complex intrapulse modulation schemes 
in a single compact formalism, we have provided a simple, physically intuitive means 
to evaluate sequence performance and inform the development of new dynamical error 
suppression protocols. 

Our approach permits an experimentalist or system designer to analytically study 
the efficacy and error-susceptibility of customized control protocols involving complex 
temporal modulation schemes (e.g. control-field phase and amplitude). For instance, 
it is possible to accurately account for error accumulation due to control imperfections 
arising from realistic constraints such as pulse overshoot, bandwidth limitations, or sim- 
ple noise in the control field itself. Further, for quantum information one may consider 
applying this approach at an algorithmic level, permitting "echo" type effects to be ex- 
ploited in minimizing the error not just of a single operation, but of a complex chain in a 
computation. With validation for the utility of our method from our study of dynamical 
decoupling, we believe that this approach will prove invaluable for future experiments 
in which it is vital to accurately and efficiently estimate the effects of real experimental 
noise on the achievable fidelity of control operation. 
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